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Abstract 

A two miniband model for electron transport in semiconductor superlattices that includes scat- 
tering and interminiband tunnelling is proposed. The model is formulated in terms of Wigner 
functions in a basis spanned by Pauli matrices, includes electron-electron scattering in the Hartree 
approximation and modified Bhatnagar-Gross-Krook collision terns. For strong applied fields, bal- 
ance equations for the electric field and the miniband populations are derived using a Chapman- 
Enskog perturbation technique. These equations are then solved numerically for a dc voltage 
biased superlattice. Results include self-sustained current oscillations due to repeated nucleation 
of electric field pulses at the injecting contact region and their motion towards the collector. Nu- 
merical reconstruction of the Wigner functions shows that the miniband with higher energy is 
empty during most of the oscillation period: it becomes populated only when the local electric 
field (corresponding to the passing pulse) is sufficiently large to trigger resonant tunneling. 

PACS numbers: 85.35.Be, 73.63.Hs, 73.63.-b, 72.20.Ht 
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I. INTRODUCTION 



Consider a n-doped semiconductor superlattice (SL) under a sufficiently large vertical 
voltage bias so that electron transport is due to resonant tunneling between minibands. For 
small voltage values, electron transport chiefly involves the lowest miniband and there are 
many appropriate kinetic theory descriptions: semiclassical Boltzmann-type equations^-, 
density matrix formulations 6,7 , transport equations for the nonequilibrium Green function 
(NGF) 8 , and Wigner-Poisson (WP) equations^. Semiclassical equations are easier to handle 
and, in particular, can be used to describe space-charge instabilities such as self-sustained 
oscillations of the current (SSOC) in dc voltage biased SLs due to the formation and dynam- 
ics of electric field domains 10 . SSOC can be found by deriving and solving a drift-diffusion 
system from the semiclassical kinetic equation^, or by a direct numerical solution of the 
latter—. Quantum transport description based on NGFs are still limited to spatially ho- 
mogeneous electric fields and therefore cannot be used to describe properly space-charge 
phenomenal. WP equations can be used to derive nonlocal drift-diffusion systems exhibit- 
ing SSOC provided collision terms are of Bhatnagar-Gross-Krook (BGK) type -. 

In contrast to work in one- miniband SL, much less is known about first-principles space- 
charge transport involving resonant tunneling in SL3. Most of the work on resonant tunnel- 
ing SL assume a large separation between time scales such that electron density and electric 
field can be assumed to be constant in each SL period and the tunneling current across 
barriers can be assumed to be stationary. Then expressions for the stationary current in an 
infinitely long SL under a constant electric field can be calculated by any quantum kinetic 
method and inserted in discrete balance equations 10 . The resulting models have been vastly 
useful to understand nonlinear electron transport in SL but they have not been derived 
from first principles. Recently, we have found a consistent perturbation method to derive 
nonlocal drift-diffusion systems (NDDS) from WP descriptions of two-miniband SLs with 
Rashba spin-orbit interaction^. However, coupling between minibands in that work does 
not contemplate resonant tunneling between them for the underlying physical description of 
the SL is too simple. 

Some time ago, Morandi and Modugno studied a variant of the standard k-p theory in 
which interband coupling terms depend on the applied electric field and used it to study 
wave function dynamics of a resonant tunneling diode 13 . For the same system, multiband 
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Wigner function approaches have also been considered^ - —. Unlu et al^ use a nonequilibrium 
Green function formulation that includes scattering due to weak coupling to a phonon bath 
to derive equations for the multiband Wigner functions. A treatment of space-dependent 
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but time-independent NGF and Wigner functions in MOSFET can be found in Ref. 
The other works focused their attention in coherent transport under an external field and 
near the semiclassical limit, thereby ignoring scattering^"—. In this paper, we present a 
simplified model of a two-miniband SL using a field dependent coupling between minibands 
similar to that introduced for resonant tunneling diodes^. We consider the corresponding 
WP system with BGK collision terms that include collision broadening and decay between 
minibands due to scattering. Electron-electron scattering is treated in the Hartree approx- 
imation through the Poisson equation. We are interested in the hyperbolic limit in which 
electric field effects, including field-dependent inter-miniband transitions, are as strong as 
the BGK collision terms and dominate electron transport. By using the Chapman- Enskog 
perturbation method, we derive nonlocal balance equations for the electron population of 
the minibands and the electric field that inherit the nonlocality of the quantum Wigner 
equation. Numerical solutions of these nonlocal equations allow us to reconstruct the time- 
resolved Wigner matrix and they exhibit resonant tunneling between minibands and SSOC. 
During SSCO, we show that the miniband with higher-energy is practically empty except 
when the local electric field is sufficiently large to allow resonant tunneling from the mini- 
band with lowest energy. Our calculations provide a first-principles description of SSCO in 
a resonant tunneling SL under dc voltage bias. 

The rest of the paper is organized as follows. Section [TT] contains the Hamiltonian we 
use as the basis of our kinetic theory. The governing WPBGK equations for the Wigner 
functions are introduced in Section IIHI The derivation of nonlocal balance equations by 
the Chapman- Enskog method is given in Section IIVI Section [V] presents numerical results 
obtained by solving the nonlocal balance equations with appropriate boundary conditions for 
the contact regions and dc voltage bias. In particular, these solutions include SSCO. Finally 
Section ED contains our conclusions and the Appendix is devoted to technical matters. 
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II. MODEL HAMILTON! AN 



Let us assume that the total Hamiltonian describing our system is 

H to tal = H + H sc , (1) 

where H sc represents scattering and H(x, —id/dx) is a 2 x 2 Hamiltonian H corresponding 
to a SL with two minibands of widths Ai and A 2 , gap energy 2g and SL period I, 

f -^-(1-coskl) -eW(x)+g eFld \ 

H(x, k) = 2 , 2) 

I eFIS Ai (1 - cos kl) -eW{x)-g j 

= [h (k) - eW(x)}(7 + h(k) ■ a + eFIS a 1} 

Here we have considered tight-binding dispersion relations for the minibands and — e < 0, 
W and — F = —dW/dx are the electron charge, the electric potential, and the electric 
field, respectively. The electric potential W in H describes electron-electron interaction in 
a self-consistent Hartree approximation. 

The matrix Hamiltonian H can be written as a linear combination of the Pauli matrices 

-A ( l 

0-0=1 I , CTl = I I , cr 2 = I , 0"3 





i / V -1 



with coefficients: 



h (k) = —a (1 — coskl), h\{k) = 0, 

h 2 {k) = 0, h 3 {k) = -7(1- cosW) +g, (3) 

_ A2-A1 _ A 2 +A! 

"45 I 4 

The term eFl5 <J\ in is a field-dependent tunneling term derived by means of the k- 
p theory for the evolution of the Wannier envelope functions [cf. Equations (33) of Ref. 
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without second order terms, i.e. with M nn > = 0]. The dimensionless parameter 5 is a 



phenomenological parameter proportional to the interminiband momentum matrix element: 

8 = -, P = - ut-—dx, (4) 

2m* gV I J_ l/2 2 dx y 1 

where 1*1,2 are the periodic parts of the miniband Bloch functions. A related model has been 
used to describe coherent transport in a resonant interband tunnelling diode^^£. 
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The miniband energies S^{k) are the eigenvalues of the free Hamiltonian Ho(k) = 
h (k)(To + h(k) ■ <x (zero electric potential), given by 

£ ± (k) = h (k)±h 3 (k). (5) 



The corresponding spectral projections are 



(6) 



so that we can write 



Ho(Jfe) = £ + (k)P + + £-(k)P-. 



(7) 



III. WIGNER FUNCTION DESCRIPTION 



If il) a (x, y, z,t), a — 1, 2, are the second quantized wave function amplitudes expressed in 
the Bloch basis, the Wigner matrix is^ 



21 r 

f ab (x,k,t) = - V / (^(x + jl/2,y,z,t)Mx-jl/2,y,z,t))e^dxj 



(8) 



where 5* is the SL cross section. Note that the Wigner matrix is periodic in k with period 
2n/l. It is convenient to write the Wigner matrix f(x, k,t) in terms of the Pauli matrices: 
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(9) 



i=0 



The Wigner components f l (x,k,t) are real and can be related to the coefficients of the 
Hermitian Wigner matrix by 



fn = f° + f, fi2 = f 1 ~if\ 

f21 = f+lf\ /22 = /°-/ 3 . 



(10) 



Hereinafter we shall use the equivalent notations 



/ 



ff°\ 

f 1 
P 



(11) 
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The populations of the minibands with energies £ ± are given by the moments: 



n (x,t) 



2tt 



it/I 



[f(x,k,t)±f(x,k,t)] dk, 



(12) 



and the total electron density is n + + n~ . 

After some algebra, from the time-dependent Schrodinger equations for wave functions ip a 
with the Hamiltonian H tot in (Q]), we can obtain the following Wigner-Poisson-Bhatnagar- 
Gross-Krook (WPBGK) equations for the Wigner components 



dp 


a 


dt 


~ h 


df 


a 


dt 





1 ■ 

— S1I 

h 



7 . 
— sn 
n 



(13) 
(14) 



whose right hand sides contain collision terms Q[f] arising from H sc . These terms will be 
modeled phenomenologically and described later. Electron-electron collisions are treated 
in the Hartree approximation and described by the Poisson equation for the electrostatic 
potential: 



d 2 W 



(n + + n - N 



D) 



dx 2 I 

where e and iV/j are the SL permittivity and the 2D doping density, respectively. In ([13 

(HI, 

2(5 — 7) + 7 cos kl A 



(15) 



h 



z7, z7= (0,0,1) 



el 



J = -CO 

00 



elS 



2 / m (x l fc,t) = -— e ijkl f?(x,t)AjF(x,t), 



el5 



J = -CO 

00 



B a f n (x,k,t) = — e ijkl fp(x,t)A+F(x,t), 



j=-oo 



©m = 0i/+ 



e 3 / 3 



(16) 
(17) 

(18) 

(19) 

(20) 



We have defined the operators 



(Afu)(x, k) = u ( x + ^,kj ± u (x — 7-, k 



(21) 



(the subscript is omitted for j = 1) and the spatial averages: 

M"J/2 



(F(x,t))j = — / F(x + s,t)ds 

J* J-jl/2 



dW 
dx 



x, t) 







A~W(x,t) 



(22) 
(23) 



dx x " " //:? jl 

12 and it contains two terms: a BGK term 



Our collision model is similar to that used in Ref. 
which tries to send the miniband Wigner function to its local equilibrium and a scattering 
term that sends electrons from the miniband with higher energy (whose electron density is 
n + ) to the miniband with lower energy (whose electron density is n~): 

f°-Q° 



Q°[f] 
Qlf] = 



f-n uf° + f 



T 



+ 



n 



v. 



m*k B T 



71 



71 



r 4 + [E - £ ± (k)Y 



fi -E 



In 1 + e k B? 



dE, 



n 



— / ^{k-n^dk. 
2tt J-k/i 



(24) 
(25) 
(26) 
(27) 

(28) 



The chemical potentials of the minibands, fi + and fi~ are calculated in terms of n + and 
rT respectively, by inserting fl271) in (|28p and solving the resulting equations. The local 
equilibria (j) are the integrals of collision-broadened 3D Fermi-Dirac distributions over the 
lateral components of the wave vector on the plane perpendicular to the growth direction x.— 
As the broadening energy r — > 0, the line-width function in the integrand of ( !27l) becomes 
5{E-S ± {k)). 

Our collision model should enforce charge continuity. To check this, we first calculate the 
time derivative of n ± using (fT2|) to f|T4l) : 



~dt 



alA- 
27xh 



2nh 



±- 



/ 

"27 



W / 1 rylA 

smkl(f°±f)dk 

tt/I 

n/l j 

e 3 fdk ' 

-tt/I 



w/l 



tt/I 



smkl(f ± f)dk 



2tx 



n 



n/l (Q°[f}±Q 3 [f})dk = T 



(29) 



where we have employed j Qif°dk = j &2f 1 dk = 0. Then we obtain: 



d_ 
di 



n 



n 



A" 



I 



■k/I 



w/l 



sin kl (af + 7 / 3 ) dk 



0. 



(30) 
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Noting that A~u(x) = I d(u(x))i/dx, we see that (TSUI) is the charge continuity equation. 
Differentiating in time the Poisson equation (j!5p . using ( 130]) in the result and integrating 
with respect to x, we get the following nonlocal Ampere's law for the balance of current: 
dF I el f w/l 

-tt/I 

Here the space independent function J(t) is the total current density. Since the Wigner 
components are real, we can rewrite (|3T|) in the following equivalent form: 



J smkl(af° + 1 f)dk) = J{t). (31) 



£ ^r + t ( a Im ^° + ^ Im fl 1 = At). (32) 



dt h 

We are using the notation /j™ for the Fourier coefficients of f m : 

oo 

f m (x,k,t) = £ fTMe m - (33) 

j=-oo 

IV. THE CHAPMAN-ENSKOG METHOD AND BALANCE EQUATIONS 

In this Section, we shall derive the reduced balance equations for our two-miniband SL 
using the Chapman-Enskog method. Note that if we were to know the Wigner matrix as a 
function of n ± and the electric field, Equations ( 129]) and the Poisson equation ( TT5l) would be 
the sought balance equations and could be solved directly. As they are now, Equations ( 129]) 
are not closed. However, in a limit in which collisions and electric potential terms dominate 
all others in the Wigner equations, it is possible to use perturbation theory to close ([29]) . 
The idea is that in this so-called hyperbolic limit, the Wigner matrix is very close to a local 
equilibrium (modified by the electric field) which depends on n and F. Using two terms 
in a Chapman-Enskog expansion, we show below that Equations ( 1291) can be closed. 

First of all, we should decide the order of magnitude of the terms in the WPBGK 
equations f lT3|) and f lT4|) in the hyperbolic limit. In this limit, the collision frequency 
1/r and the Bloch frequency eFMl/fr are of the same order, say about 10 THz. Then 
Fm = 0(h/(elr)). Typically, 2g/h is of the same order, so that the term containing 
2g/h in ( IT4l) should also balance the BGK collision term. The other terms are of or- 
der ^l/(hxo), where Xq is the characteristic length over which the field varies, and they 
are much smaller, so that A = / ~frl/(hxo) 1. From the Poisson equation, we obtain 
Xq/1 = eFmI '(eNu) = eh/(e 2 TlNi)), and therefore the small dimensionless parameter is 
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The scattering time r sc is much longer than the collision time r, and we shall consider 
t/tsc = 0(A) <C 1. Equations (|T3|) and ffT4"]) can be written as the scaled WPBGK equations 
as follows: 



L/-n = -A [r^ + A/j. 



(35) 



where we have inserted the book-keeping parameter A which is set equal to 1 at the end 
of our calculations.- 1 ^ This trick saves us from rewriting our equations in nondimensional 
units. Here the operators L and A are defined by 

ff l \ ( \ / \ 



/ 3 
\-PJ 



hf = f-TQ 1 f-TQ 2 



f° 


V o / 



re 3 



-P 
f l 
\ o / 



(36) 



A/ = V2 
where 







— sin fcZ A 



ocf + 7 



/ 3 



+ ?(cos -2) 
ft 





z7 x / 



2or r 

?7l = ~Z-, ^2 = — 



(37) 



To derive the reduced balance equations, we use the following Chapman-Enskog ansatz: 



f(x, k, t; e) = p°\k- n + , n~, F) + ^ f^(k; n + , n~, F) A m , 



m=l 



dF 



dn 
~dt 



m=0 

_|_ OO 



J2Ai(n + } n~,F) X m . 



(38) 
(39) 
(40) 



m=0 



The functions and J m are related through the Poisson equation (TT5]) . so that 



4+ _l 4- 



/ <9J„ 



e <9x 



Inserting (158" ]) to fj4"D]l into ([35]), we get 

L/« = 
L/( 2 ) = 



r 



T ■ 



dt 



A/ 



(0) 



dt 



A/ 



(i) 



a/(o) 



dt 



(41) 

(42) 
(43) 

(44) 



and so on. The subscripts and 1 in the right hand side of these equations mean that we 
replace edF/dt\ m = J6o m — J m , dn^ /dt\ m = A^, provided Sqo = 1 and 6o m = if m ^ 0. 
Moreover, inserting (I38p into (TT2"|) yields the following compatibility conditions: 



,(1)0 _ ,(1)3 

Jo ~ Jo 



(2)0 _ ,(2)3 







Jo 



(45) 
(46) 



etc. 



To solve (142]) for f^ = (p,we first note that 

oo 

-TQ X ip = i J2 '^irj''' 1 '''- 



]=~oo 



]=-oo 

CO 



-tQ 3 <p = -5 WWA+F, 



j=-oo 



rel 



Then (S2D and ([26]) yield 



Vid 2 Zj (Aj J 7 )' 



+ i- 



i 



■ V1 5 2 Z 3 (AJT) (A+F), 



-^5(1 + i^) Z d [(</>+ + 0t) zM/ATj- + 



1 + 21?- 



-?^%Mr(A+^) 



(47) 

(48) 

(49) 
(50) 

(51) 

(52) 
(53) 

(54) 



+ i- 



+ 



-rj^Z^AjT) (A+F)- 



Here we have used that the Fourier coefficients 

< r" 



COS 



7T Jo 



(Jkl) ± dk, 



are real because are even functions of k. The coefficients Z~ and M- are defined as 



Mf = — 



Z, 



1 + id* + 



8 2 {AfF) 2 
1 + id* 



(55) 

(56) 
(57) 
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The solution f(°> = <p given by (|5ip - (|54p is essentially the local equilibrium fl given by (|26p - 
(128]) modified by the field-dependent terms 9« that appear in the Wigner equations (TTBl 
and (|T4j) . This solution yields convective terms in the balance equations which contain first 
order differences. In the semiclassical limit, these equations become a hyperbolic system 
which may have discontinuous solutions (shock waves). Then it is convenient to regularize 
such solutions by keeping diffusion-like terms (second order differences) arising from the 
next-order Wigner functions 

The solution of ( 143]) is p 1 ' = ?p with 



1 + ify 



1 



S 2 M+{AjTf 



(58) 



1 



, 1 _ bA]T 
M + r + r H r 6 



77! (1 + M+M, 



77! (1 + M+M: 



M , r f + 1 -q r ? + r f + i ^7T r f 

J 1 + J J 1 • n), 3 



M~r z 

3 3 



1 + ii£ ^ ~~ ^ ~ 1 + ^ Tj 



1 + itf, 



5 2 M"(A+^) 2 



(59) 
(60) 
(61) 



5A+.F 



1 + J 



j 



J7i(l + i^-) (1 + M+Mr 
Here r is the right hand side of ( H3l) . 

The balance equations can be found in two ways. We can calculate for m — 0, 1 in 
( H0|) by using the solvability conditions (|4"5]) and (jl6"j) in (H3|) and (jHJ), respectively. More 
simply, we can obtain the balance equations by inserting the solutions fl5Tl) to (]54"P and fl58|) 
to (I6T]) in the balance equations ( 129]) and in the Ampere's law (1311) . The result is: 



dt 

at 



n 



F) = T R(n + ,n-,F), 



e — + e (D + (n + , n~,F) + D.(n + , n~ , F)) x = J(t) 
a ± 7 



± 



Im(^±^ + ^°±^i), 



= - [r] 2 n + + 28T{<pl + <A 



(62) 
(63) 
(64) 
(65) 



Note that Eq. (163]) can be obtained from (162]) and the Poisson equation. Equations (162]) to 
(165]) must be solved together with the Poisson equation (fT5]) . the expression for the local 
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equilibrium Wigner densities ( |2"T|) and expressions ( 128|) for n . The zeroth and first order 
Wigner functions cpj and ipj in and f )65|) can be obtained from Equations fl5T]) - fl5^|) and 
IBTi). respectively. The complete expressions for £)± and i? are shown in Appendix fAI 



V. NUMERICAL RESULTS 

To solve numerically the system of equations ( 1621) - ( )65l) . we have to add the voltage bias 
conditions for the electric potential and appropriate boundary conditions at the contact 
regions. Note that our equations involve finite differences and several one-period integral 
averages. This means that we need to give boundary conditions over intervals of size 21 
before x = and after x = Nl, not just boundary conditions at x = 0, Nl as we would 
give for semiclassical drift-diffusion equations. At the injecting region (cathode), the usual 
boundary condition is that the electron current density satisfies Ohm's law and therefore it 
is proportional to the electric field there. We use this condition for each point of the interval 
—21 < x < 0. Similarly, we also need the electron densities n at the cathode. To avoid 
inconvenient boundary layer effects, we choose their values for a spatially uniform stationary 
state with a given value of the field. The resulting boundary conditions in —21 < x < are: 
W = and 

dF 

£ ~fo + ^cathode F = J, (66) 

n ± = n± (67) 

where nf t are the miniband electron densities corresponding to a spatially uniform stationary 
state. The latter can be obtained by equating to zero the right hand sides of the rate equation 
(1621) and the Poisson equation ( TT5l) : R(n + , n~ , F) = and n + + n~ = Nd, respectively. The 
result is 

2 t2\ 



712 ~ Nd {2 T 85^ + 2 V2 (1 + V f + WW) J ■ {Q8) 
The boundary conditions in the anode region (Nl < x < Nl + 21) are: W = V and 

dF ,n + + n~ . , . 

S -Qjr + °*node ( ) F = J, (69) 

n + = 0. (70) 

The lower miniband electron density rT in the anode region is obtained from the Poisson 
equation f lT5|) . 
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FIG. 1: (a) Electron current vs field in a spatially uniform stationary state, (b) Total current 
density vs time, (c) Electric field profile at different times of one current self-oscillations cycle. 
At time (1) the field is above the resonant value for the middle SL point x = Nl/2. (d) Electron 
densities n ± /(n + + n~) vs time for point x = Nl/2. When the electric field is above the resonant 
value (time (1)), the electron transport between minibands occurs. 



To present numerical results, we have used the parameter values corresponding to a 
GaAs/AlAs SL from Table I of |l9| which has narrow minibands so that resonant tunneling 
plays an important role in electron transport. Our parameter values are: ds = 1.5 nm, 
d w = 9 nm, I = d B + d w = 10.5 nm, N D = 2.5 x 10 10 cm" 2 , r = 0.0556 ps, r sc = 0.556 
ps ,20 v = 9 V) N = 200, (y cathode = 1.4 fi^m" 1 , a anode = 0.7 fi^nr 1 , T = 5K,A 1 = 2.6 
meV, A 2 = 13.2 meV, P/h = 0.2238/nm,^- T = 1 meV.^ With these values, a = 2.6 
meV, 7 = 3.9 meV, 5 = 0.12. We have selected the following units to present our results 
graphically: F M = hj{eW) = 11.28 kV/cm, x = eF M l/(eN D ) = 31.4 nm, t Q = h/a = 0.25 
ps, J = aeN D /h = 1.58 x 10 4 A/cm 2 . 

Figure [H (b) illustrates the resulting stable self-sustained current oscillations. They are due 
to the periodic formation of a pulse of the electric field at the cathode x = and its motion 
through the SL. Figure Q] (a) depicts the electron current vs field in a spatially uniform 
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-1 -0.5 0.5 1 -1 -0.5 0.5 

kl/% kl/n 



FIG. 2: (a)-(b) Wigner matrix off-diagonal terms f 1 and f 2 vs k, at time (1) (tunneling transport), 
and time (2) (no tunneling), (c)-(d) Wigner matrix diagonal terms / ± f 3 vs k 

stationary state, with a local maximum at the field resonant value 2g/(el). Figured] (c) 
depicts the electric field profile at different times during one self-sustained current cycle. 
Figure [T] (d) shows the tunneling transport between minibands when the electric field is 
above the resonant value (time (1)) calculated at the middle point of the SL (x = Nl/2). 

Figure [2] shows the Wigner matrix elements f \ from equations (|5TT)-()54l). (!58l) -( l6TT) and 
( |33|) . for the middle SL point (x = Nl/2) vs k at times (1) (with tunneling transport 
between minibands) and (2) (with no tunneling). Figure |2fa)-(b) illustrates the Wigner 
matrix off-diagonal terms f 1 and f 2 , which are responsible for the tunneling transport 
between minibands. Figure E)^c)-(d) shows / ± / 3 , which are related with the electron 
densities n^. 

Figure |3] illustrates the effect of varying the voltage bias on the total current for a iV = 60 
period SL. Figure |3] (a) depicts the total current density average, maximum and minimum 
values for different voltages. It can be seen that when the bias is above a critical voltage, 
the current self sustained oscillations appear and their amplitude increases from zero at the 
bifurcation point. This circumstance does not depend on whether the voltage is increasing 
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FIG. 3: (a) Total current density (average, maximum and minimum values) vs voltage bias, (b) 
Current oscillation frequencies vs voltage bias. 



or decreasing, therefore the critical voltage corresponds to a supercritical Hopf bifurcation. 
Figure |3] (b) shows that the oscillation frequencies decrease as the voltage increases above 
its critical value. Immediately above the critical voltage, self-oscillations are due to repeated 
triggering of small pulses of the electric field that die near the cathode and before they can 
reach the end of the SL. As the voltage increases, the pulses are able to grow and reach the 
anode region. Since their average velocity does not vary that much, the oscillation frequency 
is correspondingly smaller. In a transition region between 1.5 and 3V, the current oscillation 
is somewhat irregular. The region of self-oscillations ends at a larger voltage of about 5.3V. 
Similar phenomena are observed in models of the Gunn effect in bulk GaAs. See Chapter 6 



in Ref. 
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If we use parameters corresponding to a weakly coupled SL with miniband widths below 
1 meV (that come from using wider quantum barriers), we run into problems of numerical 
convergence and, possibly, breakdown of the Chapman- Enskog perturbation scheme. To 
explore the limit of weakly coupled SL, a different perturbation scheme based on miniband 
smallness seems necessary. This is outside the scope of the present paper. 



VI. CONCLUSIONS 



For strongly coupled SLs having two populated minibands, we have introduced a kp 
Hamiltonian that contains a field-dependent tunneling term and derived the corresponding 
Wigner-Poisson-BGK system of equations. The collision model comprises two terms, a BGK 
term trying to bring the Wigner matrix closer to a broadened Fermi-Dirac local equilibrium 



15 



at each miniband, and a scattering term that brings down electrons from the upper to the 
lower miniband. By using the Chapman-Enskog method, we have derived quantum drift- 
diffusion equations for the miniband populations which contain generation-recombination 
terms. As it should be, the recombination terms vanish if there is no inter-miniband scatter- 
ing and the off-diagonal terms in the Hamiltonian are zero. These terms represent miniband 
coupling due to the electric field and originate the resonant tunneling transport. For a 
superlattice under dc voltage bias in the growth direction, numerical solutions of the cor- 
responding quantum drift-diffusion equations show self-sustained current oscillations due to 
periodic recycling and motion of electric field pulses, and resonant tunneling between mini- 
bands when the electric field is above the resonant value. Numerical reconstruction of the 
Wigner functions during self-oscillations confirms this picture. 
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Appendix A: Detailed expressions for D± and R 

The recombination term R(n + , n~, F) (1651) depends on (p^ and which can be obtained 



from (1531) and ( 1601) for j = 0, taking into account that A J 7 = 0, Aq J 7 = 2J-" and 



'0 



n 



^0 



n + — n 



l + r) 2 + ^ 2 T 2 
1 

1 + nj + A5 2 F 2 
^T 2 



ar 



' 5(1 +ri 2 - A5 2 F 2 )(l - ri 2 - A5 2 F 2 (n + - n~)) 37_ 
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f dn + dn~ 



dn 



dt 



l-n 2 - 45 2 F 2 (n + - 
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The time derivatives 



drr 11 



dt 



Chapman-Enskog expansion o 



and 



dF 
~dt 



are obtained from the first two terms of the 



(1521) and I^H]) respectively: 



dt 



dF 
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A-Df = =F— T 
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+ 



1 + tf 2 



, r(l + r/ 2 + 4<RF 2 ) 
+ th5 2 (TA-J r A + TReZ 1 + 



ReZ 1 ((A"J r ) 2 ImM+ + (A + J) 2 Imikff ) + 
ImZ 1 ((A-J r ) 2 ReM+ + (A + J^) 2 ReMf ))) + 
^77i5 2 [ReZi((A _ J r ) 2 ImM 1 + - (A+Jj^mMf 
ImZ 1 ((A-J r ) 2 ReM+ - (A+F) 2 ReMf)]] 
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The expression of _D±(n H 
expansion D± and d£ . 



n 



F) is based on the first two terms of the Chapman-Enskog 



D±(n + ,n-,F) = D±(n + , n~, F) + D±(n + , n~, F) 
Where D±^ and follows: 

D ( £\n\n-,F) = ^bn(v? ± = 

(t^| T^ 1 5 2 (A-^)(A + ^)ReZ 1 

+ ^!((I m M+(A-J^) 2 + ImMf (A+F) 2 )ReZi 
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The functionals S 1 (n + ,n-,F), Af(F), Bf(F) and Cf(F) are as follows: 
1 
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Where the real and imaginary parts of ipj are: 
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- (f)j)A+F(ReM-(ImZ j + VjReZj) + lmM~{ReZ j - fylmZj))] 

2 - 7]i5 2 (A+J r ) 2 (ReM7ReZ i - ImAfrimZ,-) 



Re<^ 3 = 



1 



1 + tfj 



2 



2 

^ 2 



1 + tf 2 



?7i5 2 (A+J r ) 2 (ReM-ImZ i + ImAfrReZ,-) 



+ 



(0+ + 0-)A+^A7^ReZ j 



19 



Now we can obtain the expressions for 
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, dn + 



dt 



dn 



dt 



x 



^(A+J 7 ) 2 (ReMfReZi - ImMflmZi) 



+ 



-2i?i 



2^5 2 A+^A+ ^ 



(1+^)2 ft 

(ReMfReZi - ImMflmZi) 



*7i<f 2 (A+F) 2 

aimZi 



/ 9ReMf 



V at 



ReZi - ReMf 



dReZi 



at 



\m.Zi 



dImM{ 



dt 



lmM[ 



dt 



o 

, 9n 



(tf + *r)(A+ — 

aimZi 



at 

a^ 



at 



A+.FA "JlmZi+ 



at 



at 



dT 

A~J r lmZ 1 + A- — 
dt 



A+J^ImZi + 



o 

, dn + 



dt 



<9n 



at 



x 



^(A+J 7 ) 2 (ReMflmZi + ImMfReZi) 



+ 



$\ - 1 &&! 



2r ll 5 2 A + FA 



(1+tf 2 ) 2 ^ 
a^ 



mS 2 (A + Tf 



dt 

( dReM{ 



(ReMflmZi + ImMfRe^) - 



ImZi + ReMf 



aimZx 



at 



+ ReZ 1 



aimMf 



at 



+ 



ImMj" 



_ aReZi 



at 

(<t>t: 





, dn + 



+ 



dt 



+ (0D 



, dn' 



dt 



A^J'ReZi + A 



) A+J r A-J r ReZ 1 + 

0/ 



at 



A+J^ReZi + 



aReZi 



at 



In the above expressions we have used (0f)' = d^/dn^ and dA/dt\ = {dJ r /dt\ ) v We 
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also need to calculate Zj, Mp, dZ 1 /dt\ and dM^/dt\ 
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The functionals Mf{F) are as follows: 
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Finally, we need to calculate the time derivative 



dt 
OImM± 



dt 



dt 



25 2 A ± J r 

1 

Vi 



OA 



dt 



A ± — 

dt 



(1+tf 2 ) 2 



^iA ± J r dd l 



l + dt 
1 - ■&j)A ± J : 



ddx 



dt 



dT 

+ 20 1 (l + *i)A ± ^ 



1 S.A. Ktitorov, G.S. Simin and V.Ya. Sindalovskii, Sov. Phys. Solid State 13, 1872 (1972). 

2 A. A. Ignatov and V.I. Shashkin, Sov. Phys.-JETP 66, 526 (1987). 

3 A. A. Ignatov, E.P. Dodin and V.I. Shashkin, Mod. Phys. Lett. B 5, 1087 (1991). 

4 A. Sibille, in Semiconductor superlattices: growth and electronic properties, ed. by H.T. Grahn. 
Page 29 (World Sci., Singapore, 1995). 

5 L.L. Bonilla, R. Escobedo and A. Perales, Phys. Rev. B 68, 241304(R) (2003). 

6 V. V. Bryksin and P. Kleinert, J. Phys.: Cond. Matter 9, 7403 (1997). 

7 M. V. Fischetti, J. Appl. Phys. 83, 270 (1998). 

8 H. Haug and A.-P. Jauho, Quantum Kinetics in Transport and Optics of Semiconductors. 2nd 
ed (Springer, Berlin 2008). 

9 L.L. Bonilla and R. Escobedo, Math. Mod. Meth. Appl. Sci. (M 3 AS) 15, 1253 (2005). 

L.L. Bonilla, H.T. Grahn, Rep. Prog. Phys. 68, 577 (2005). 

1 E. Cebrian, L.L. Bonilla and A. Carpio, J. Comput. Phys. 228, 7689 (2009). 

2 L.L. Bonilla, L. Barletti and M. Alvaro, SIAM J. Appl. Math. 69, 494 (2008). 

3 O. Morandi and M. Modugno, Phys. Rev. B 71, 235331 (2005). 

4 M.B. Unlu, B. Rosen, H.L. Cui and P. Zhao, Phys. Lett. A 327, 240 (2004). 

5 L. Demeio, L. Barletti, A. Bertoni, P. Bordone, and C. Jacoboni, Physica B 314, 104 (2002). 

6 O. Morandi and L. Demeio, Transport Theor. Stat. Phys. 37, 437 (2008). 

7 O. Morandi, Phys. Rev. B 80, 024301 (2009). 

8 H. Jiang and W. Cai, J. Comput. Phys. 229, 4461 (2010). 

9 J. Kastrup, R. Hey, K. Ploog, H.T. Grahn, L. L. Bonilla, M. Kindelan, M. Moscoso, A. Wacker 
and J. Galan, Phys. Rev. B 55, 2476 (1997). 

;0 J. Shah, in Hot Carriers in Semiconductor Nanostructures: Physics and Applications, pp. 279 
(Academic Press, Boston 1992). 

25 



For numerical estimates, see also W. T. Wenckebach, Essential of Semiconductor Physics (J. 
Wiley & Sons, Chichester, 1999). 

L. L. Bonilla and S. W. Teitsworth, Nonlinear wave methods for charge transport (Wiley- VCH, 
Weinheim, 2010). 



26 



